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Abstract —  This  article  presents  one  modification  of 
algorithm  for  generation  of  n-dimensional  rotation  matrix 
M,  which  rotates  given  n-dimensional  vector  X  to  the 
direction  of  given  n-dimensional  vector  Y.  Algorithm, 
named  N-dimensional  Rotation  Matrix  Generation 
Algorithm  (NRMG)  includes  rotations  of  given  vectors  X 
and  Y  to  the  direction  of  coordinate  axis  Xi  using  two- 
dimensional  rotations  in  coordinate  planes.  Proposed 
modification  decreases  the  number  of  needed  two- 
dimensional  rotations  to  2(Lw-l)  were  Lw  is  the  number 
of  corresponding  components  of  the  two  given  vectors, 
that  are  not  equal. 
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I.  INTRODUCTION 

High-dimensional  spaces  frequently  occur  in  mathematics 
and  the  sciences,  in  example  n-dimensional  feature  space, 
which  presents  input  signals  of  neural  network  or 
collection  of  n-dimensional  parameters  for 
multidimensional  data  analysis.  Rotation  is  one  of  rigid 
transformations  in  geometrical  space,  which  preserves 
length  of  vectors  and  can  be  presented  using  matrix 
operation  like  Y  =  M.X  ,  where  X  and  Y  are  input  and 
output  vector  respectively,  and  M  is  rotation  matrix.  This 
article  proposes  one  modification  of  n-dimensional 
rotation  matrix  generation  algorithm,  which  increases 
performance  if  some  of  corresponding  components  of  the 
two  given  vectors  are  equal. 

H.  SHORT  DESCRIPTION  OF  ALGORITHMS 
FOR  N-DIMENSIONAL  ROTATION  MATRIX 
GENERATION 

Let's  say  that  we  have  two  n-dimensional  vectors  X  and 
Y,  having  the  same  dimension,  X,  Y  e  R".  We  want  to 
obtain  a  rotation  matrix  M  that  satisfies  the  equation. 

Y  =  M.X  (1) 

where  Y  has  the  same  norm  as  X  and  the  same  direction 
as  Y,  i.e.  ||y||  =  ||X||  and  cos(y,  y)=  1 . 


One  of  possibilities  to  generate  rotation  matrix  M  is  to 
rotate  given  vectors  in  two-dimensional  subspace, 
generated  by  vectors  and  map  it  back  to  Rn  as  follows 
[4], [19]: 

Obtain  two  orthogonal  vectors  u=X/||X||  and 
v  =  (Y  —  (u.y)u)/||(Y  —  (u.y)u)||  that  have  the  same 

norm.  Then  P=uuT+vvT  is  a  projection  onto  the  space, 
generated  by  X  and  Y,  and  Q=I  -  uuT  -  vvT  is  a  projection 
onto  the  n-2  dimensional  complementment  subspace.  So 
the  rotation  just  has  to  take  place  on  the  range  of  P.  This 
rotation  has  to  be  performed  in  R2.  Then  result  is  mapped 
back  to  Rn  by  [a,b]T— >au+bv.  So  matrix  of  rotation,  which 
rotates  X  to  the  direction  of  Y  is  calculated  as  follows: 


M  =  I  -  uu  T  -  wT  +  [u  v]R(e)[u  v]T  (2) 


where  R(0): 


cos 


(O)  -  sin  (9) 


and  cos(0)  = 


XTY 


sin(0)  cos(9)  v  7  ||X||||Y|| 

Another  way  to  generate  rotation  matrix  M  is  to  use 
Householder  Reflection  [5],  [11],  [18],  [19].  If  u  and  v  are 
vectors  with  the  same  norm,  in  example  obtained  from 
input  vectors  X  and  Y  as  u=X/||X||  and  v=Y/||Y||, 
there  exists  an  orthogonal  symmetric  matrix  P  such  that 


v  =  P.u 


where 


P  =  I  -  2W.W J 


and 


W  =  (u  —  v)/||u  —  v|| .  Matrix  P  is  matrix  of  reflection 

(not  a  rotation)  because  det  P  =  -1,  (which  gives  the  name 
to  method)  that’s  why  to  obtain  matrix  of  rotation  M,  for 
which  det  M  =  1,  have  to  be  performed  two  subsequent 
reflections.  Matrix  of  rotation  can  be  obtained  as 
multiplication  of  two  matrices  of  reflection  PI  and  P2  as 
M  =  P1.P2  . 

The  third  way  to  generate  rotation  matrix  M  is  to  use 
rotation  of  given  vectors  to  the  direction  of  one  of 
coordinate  axes  (i.e.  axis  3q ).  Algorithm,  named  N- 


dimensional  Rotation  Matrix  Generation  Algorithm 
(NRMG  algorithm)  [1]  includes  the  following  sequence 
of  operations: 

1)  Obtaining  rotation  matrix  Mx,  which  rotates  given 
vector  X  to  the  direction  of  axis  x\  as  follows: 
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n-1 

Mx  =  nG(n-k’n-k  +  1’0n-k)  (3) 

k=l 


where  coefficients  sin(0k)  and  cos(0k)  of  Givens 
matrices  G(k,k  +  1  ,  0k)  are  calculated  using  formulas: 


sin(0k)  =- 


cos(0k)  = 


xk+l 


xk  +  xk+l 


xk 


xk  +  xk+l 


tfxk  +  xk=l>0 


(4) 


sin(  0k)  =  0,  cos(0k)  =  1  if  xk  +  xk=i  =  0 


2)  Obtaining  rotation  matrix  My,  which  rotates  given 
vector  Y  to  the  direction  of  axis  using  (3)  and  (4). 

3)  Obtaining  rotation  matrix  M  as  multiplication  of  Mx 
by  inverse  matrix  of  My  given  as  M=Mx  '  -My  [20], 


III.  DEFINITION  OF  THE  TASK 

Let's  say  that  we  have  two  n-dimensional  vectors  X  and 
Y,  having  the  same  dimension  and  the  same  norm,  X,  Y 
g  Rn  ,  |Y|  =  ||X||  and  let’s  some  of  their  corresponding 

components  Xj  and  y,  are  equal:  x;=y;  if  ieK<z{l,...n}. 
We  want  to  obtain  a  rotation  matrix  M  that  satisfies  the 
equation. 

Y  =  M.X  (5) 

As  it  can  be  seen  from  this  definition,  proposed 
modification  of  NRMG  algorithm  needs  the  two  given 
vectors  X  and  Y  to  have  the  same  norm,  while  the 
original  NRMG  algorithm[20]  do  not. 


IV.  DESCRIPTION  OF  MODIFIED  NRMG 
ALGORITHM 

Modified  NRMG  (MNRMG)  algorithm  for  generation  of 
N-dimensional  rotation  matrix  M  is  NRMG  algorithm  [1], 
which  uses  rotations  in  coordinate  planes  only  in 
subspace  V,  that  is  spanned  by  the  unit  vectors  of  axes,  on 
which  components  of  given  vectors  X  h  Y  are  not  equal. 
MNRMG  algorithm  contains  the  following  operations: 

1)  Comparing  corresponding  components  of  the  two 
given  vectors  X  h  Y  and  obtaining  vector  ,  which 
consist  of  indexes  of  components  that  are  not  equal. 

2)  Obtaining  rotation  matrix  Mx  ,  which  rotates 
projection  of  given  vector  X  in  subspace  V  to  the 
direction  of  one  of  axes,  having  indexes  in  w  (i.e. 
the  first  one) 

3)  Obtaining  rotation  matrix  My,  which  rotates 
projection  of  given  vector  Y  in  subspace  V  to  the 
direction  of  the  same  axis  as  for  projection  of  X. 

4)  Obtaining  rotation  matrix  M,  which  rotates  given 
vector  X  to  the  direction  of  given  vector  Y  as 


multiplication  of  matrix  Mx  and  inverse  matrix  of 
My  ,  as  follows: 

M=My  .Mx  =My.Mx  (6) 

Multiplication  of  vector  by  matrix  Mx  or  by  matrix  My  do 
not  change  vector’s  components,  which  indexes  are  not  in 
vector  of  indexes  w  .  Thus  the  length  Lw  of  the  vector  w 
defines  the  number  of  rotations  in  coordinate  planes  (base 
operations  of  algorithm,  described  below),  which 
coefficients  have  to  be  calculated  to  obtain  searched 
matrix  M.  Below  operations  of  proposed  algorithm  will 
be  described  subsequently 

4.1.  Comparing  corresponding  components  of 

THE  TWO  GIVEN  VECTORS  X  H  Y  AND  OBTAINING 

VECTOR,  WHICH  CONTAINS  INDEXES  OF  COMPONENTS 
THAT  ARE  NOT  EQUAL 

To  obtain  vector  w  ,  which  contains  indexes  of 
corresponding  components  of  the  two  given  vectors  X 
and  Y,  that  are  not  equal,  have  to  be  performed 
subsequent  comparing  of  corresponding  components  of  X 
and  Y.  Indexes  of  components,  that  are  not  equal,  are 
saved  subsequently  in  vector  w  .  In  example  if  X=[l  2  3 

4]T  and  Y=[l  3  2  4]T  then  after  comparing  of 
corresponding  components  vector  w  will  have  content 
w  =[2  3]t.  where  2  and  3  are  indexes  of  second  and  third 
components,  that  are  not  equal.  The  length  Lw  of  vector 
w  is  equal  to  the  number  of  different  corresponding 
components  of  vectors  X  h  Y. 

Vectors,  that  have  zero  values  of  components,  which 
indexes  are  not  in  w  ,  generates  subspace  VcR".  In 
example  vectors  X2=[0  2  3  0]T  and  Y2=[0  3  2  0]T, 
obtained  from  vectors  X  and  Y,  given  above,  belongs  to 
subspace  VcR4.  It  is  easy  to  see  that  X2  and  Y2  are 
projections  of  vectors  X  and  Y  onto  subspace  V. 

4.2.  Obtaining  rotation  matrices  Mx  and  My, 

WHICH  ROTATES  PROJECTIONS  OF  GIVEN  VECTORS  X 
AND  Y  IN  V  TO  THE  DIRECTION  OF  ONE  OF  AXES, 
HAVING  INDEXES  IN  (I.E.  THE  FIRST  ONE) 
Rotation  of  projection  of  given  vector  X  in  V  to  the 
direction  of  the  first  of  coordinate  axes,  having  indexes  in 
w  can  be  performed  by  subsequent  multiplications  by 
Givens  matrices  [8]  as  follow: 

XI  =  nG(wk>wk+L0wk)-x  (7) 
k=Lw-l 

where  wk,wk+jare  components  of  vector  w  .  Givens 
matrices  G(wk,  wk+1,  0wk),  k  =  Lw-1,  Lw-2,  ...  1  are 
defined  as  follows  [8],  [13]: 
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G(wk,wk+1,eWk  )= 

1  0  •••  0  •••  0  •••00 

01...  0  •••  0  •••00 


0  0 


-wk 


0  0  ■••  st 


-Swk  •••  o  0 


Cwk  0  0 


(8) 


0  0  •••  0  •••  0  •••10 

0  0  •••  0  •••  0  •••01 

where  0wk  is  the  angle  of  rotation  and  coefficients 
Cwk=cos(0wk)  and  Swk=sin(0wk)  appears  at  the 
intersections  of  wk-th  and  wk+l-th  rows  and  columns.  As 
it  is  shown  in  [1]  coefficients  Cwk=cos(0wk)  and 
Swk=sin(0wk)  can  be  calculated  using  formula: 
x, 

sin(0w. 


.)=- 


L  w 


k+1 


2  2 

x  +  x 
wk  wk+1 


cos(0w,  )  = 


xk 


x2  +x2 
wk  wk+1 


if  x”  +  x~  >0 
wk  wk+1 


sin(  0 


w.  )  =  0,  cos(0w 


,  )  =  1  if  x2  +x2  =0 

k  wk  wk+1 


(9) 

Every  one  multiplication  of  vector  by  Givens  matrix 
G(wk,  wk+1,  0wk)  performs  rotation  of  its  projection  in 
coordinate  plane  (xwk  ’xwk+j )  >  which  changes  values 

only  of  vector’s  coordinates 


xwk’xwk+i 


to 


.This  multiplication  can  be  presented  as 
by  sub  matrix 


xwk’xwk+i 


multiplication  of  coordinates 
A(wk,  wk+1, 0wk)  as  follows: 

xwk 
XWk+l 


(11) 


Aw 


k  •  wk+l  ■  ®wk 


xwk 

Xwk+1 

cos(6Wk)  -sin(0Wk) 
cos(6Wk ) 


■wk 

;wk+l 


sm(6Wk) 

Taking  in  consideration  that  multiplication  with  Givens 

matrix  G(wk,  wk+1,  0wk)  changes  only  values  of  vector’s 

coordinates  xWl  ,  x ,,, ,  ,  (to  xw,  ,xw,  ,  ),  schema 

wk  wk+l  wk  wk+l 

of  multiplication  by  Givens  matrix  G(wk,  wk+1,  0wk)  can 

be  presented  as  operator  for  two-dimensional  rotation  as 

follows: 


Xi 


x2  - 


-Xi 

-x2 


Xw  k  Xw  k 

Xwk+1  *T/6w>T^wk+l 


Xw, 


'“'"'fc+l 


■Xw, 


^"k+1 


xn-l 


-Xn-1 

-Xn 


Fig.  1  Schema  of  two-dimensional  rotation,  performed  by 
Givens  matrix  G(wk,  wk+1,  Qwk) 

The  target  of  multiplication  by  Givens  matrix 
G(wk,  wk+1,  0wk)  is  to  set  to  zero  coordinate  xwk+j  •  H 

is  easy  to  find  that  equation  xwk+j  =  0  and  equations 

(10)  are  satisfied  simultaneously  when  sin(0wk)  and 
cos(0wk)  are  calculated  using  formula  (11).  Thus 
searched  matrix  Mx  can  be  calculated  as  multiplication  of 
Givens  matrices  as  it  is  given  in  (7).  Calculation  of  angles 
of  two-dimensional  rotations  0k,  k=Lw-l,  Lw-2,...,1 
really  is  not  needed.  Schema  for  rotation  of  vector  X  to 
the  direction  of  the  first  of  axes,  which  indexes  are  in  w 
using  two-dimensional  rotations  can  be  presented  as 
follows: 


'  =  •/? 


+  xr+...4  xr 


Fig.  2  Schema  for  rotation  of  vector  X  to  the  direction  of 
first  of  axes,  having  indexes  in  w 

If  calculation  uses  parallel  execution  of  two-dimensional 
rotations  [4], [16]  (named  Accelerated  Rotation  AR  in 
[20])  matrices  Mx  and  My  are  calculated  as  multiplication 
of  T log2Lw~|  matrices  of  stages  for  which  coefficients  Cwk 
and  Swk  are  calculated  using  (9).  The  proposed 
modification  of  NRMG  algorithm  avoids  execution  of 
base  operations  for  vector’s  components,  which  indexes 
are  not  in  vector  w  .  In  example  let’s  we  have  two  8- 
dimensional  vectors  X  and  Y,  that  have  two  equal 
components,  i.e.  X3=y3  and  X6=y6,  Then  schema,  which 
perform  rotation  of  X  to  the  direction  of  Y  using 
MNRMG  algorithm  and  AR  will  look  as  follows: 
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Fig.  3:  Schema  for  rotation  of  8-dimensional  vector  X  to  the  direction  of  vector  Y  using  Modified  NRMG  algorithm  and  AR  if 

Y  =  X  ,  X3=y3  and  X6=y6 


Every  one  of  matrices  Mx  and  My  is  calculated  by 
multiplication  of  Tlog2Lwl  matrices  of  stages  (in  this 
example  Tlog26]=3)  for  which  coefficients  Cw  and  Sw  are 
calculated  using  (9).  It  is  important  to  note  that  Modified 
NRMG  algorithm  needs  that  two  given  vectrors  X  and  Y 
have  the  same  norm  |Y|  =  ||X||  while  the  original  NRMG 
algorithm  [20]  do  not. 

V.  NUMERICAL  EXPERIMENTS 

Has  been  created  Matlab  program  to  test  proposed 
Modified  NRMG  algorithm.  Program  contains  function 
for  accelerated  rotation  fnMAR  and  code,  that  uses  this 
function  to  obtain  matrix  M,  which  rotates  given  vector  X 
to  the  direction  of  given  vector  Y. 
function  R  =  fnMAR(X,w) 

N  =  length(X);  %X  have  to  be  row  vector  (transposed) 
Lw=  length(w);  %w  have  to  be  row  vector  (transposed) 
if  N<=Lw  %  Length  of  X  can't  exceed  the  length  of  w 

R=  eye(N);  %  Initial  rotation  matrix  =  Identity  matrix 
step  =  1 ;  %Initial  step 

while(step<N)  %Loop  to  create  matrices  of  stages 
A=  eye(N); 
n=l; 

while(n<=N-step  &&  w(n+step)  >  0) 
r2  =  X(w(n))*X(w(n))  + 

X(w(n+step))*X(w(n+step)); 
if  r2  >  0 
r  =  sqrt(r2); 

pcos  =  X(w(n))/r;  %Calculation  of  coefficients 


psin  =  -X(w(n+step))/r; 

%  Base  2-dimensional  rotation 
A(w(n),  w(n))  =  pcos; 

A(w(n),  w(n+step))  =  -psin; 

A(w(n+step),  w(n))  =  psin; 

A(w(n+step),  w(n+step))  =  pcos; 

X(w(n+step))=0; 

X(w(n))=r; 

end; 

n=n+2*step;  %  Move  to  the  next  base  operation 
end; 

step  =  step*2; 

R=  A*R;  %  Multiply  R  by  current  matrix  of  stage  A 
end; 
end; 

Example  1  Code  of  Matlab  function  for  accelerated 
rotation  of  vector  X  to  the  direction  of  first  axis  having 
index  in  w  .  Function  returns  matrix  of  rotation  Mx 
Matlab  code,  that  uses  this  function  to  obtain  matrix  M, 
which  rotates  given  vector  X  to  the  direction  of  given 
vector  Y  is  given  below: 

N  =  length(X);  %X  and  Y  have  to  be  row  vectors 
normX=  norm(X);  %(transposed) 
normY=  norm(Y); 

if  normX  ~=  normY  %Set  norm  of  Y  equal  to  norm 
Y=(  normX/normY)*Y;  %  of  X  if  they  are  different 
end; 

w=zeros(l,N);  m=l;  %  Initialization  of  vector  w 
for(n=l  :N)  %Loop  to  create  vector  of  indexes  w 
ifX(n)  ~=  Y(n) 
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w(m)=n;  %save  in  w  index  of  not  equal  elenents 
m=m+l; 
end 
end 

Mx=fnMAR(X,w); 

My=  fnMAR(Y,w); 

M=  My'*Mx; 

Z=M*X'; 

if  max(abs(Z’-Y))  <  10A-6 
disp('Z  and  Y  are  identical'); 
end 

Example  2  Matlab  code,  which  creates  matrix  of  rotation 
M  using  fnMAR  function 

For  test  data  have  been  used  the  following  two  images: 
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Images  can  be  presented  as  64-dimensional  vectors  as 
follows  (written  using  Matlab  Language  syntax  for  row 
vectors): 

X  =[0  000000000011000001001000010 
0100001111000010010000100100000 
0000  0]; 

Y  =[0  000000000111100001001000010 
0100001111000010000000100000000 
00000]; 

Test  program,  given  above,  obtain  vector  Z  to  the 
direction  of  vector  Y  using  Modified  NRMG  algorithm. 
Program  displays  message  “Z  and  Y  are  identical”,  which 
shows  that  vectors  Z  and  Y  are  identical  at  least  with 
precision  of  10'6.  In  this  example  Modified  NRMG 
algorithm  decreases  the  number  of  executed  base 
operations  to  4  (which  is  the  number  of  different 
corresponding  components  of  the  two  vectors  X  and  Y) 
while  for  the  same  vectors  NRMG  algorithm  without 
modification  executes  70  base  operations. 

VI.  CONCLUSION 

In  this  article  we  proposes  one  modification  of  NRMG 
algorithm  for  generation  of  n-dimensional  rotation  matrix 
M,  which  rotates  given  n-dimensional  vector  X  to  the 
direction  of  given  vector  Y  having  the  same  dimension 
and  the  same  norm.  Modification  decreases  the  number  of 
base  operations  to  2(Lw-l)  where  Lw  is  the  number  of 
corresponding  components  of  given  vectors,  that  are  not 
equal.  Thus  when  the  two  given  vectors  X  and  Y  have 
significant  number  of  corresponding  components  that  are 
equal.  Modified  NRMG  algorithm  can  be  more  effective 
than  existing  algorithms,  shortly  described  above. 

It  is  important  to  note,  that  no  one  of  algorithms  for  n- 
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dimensional  matrix  generation,  described  in  p.2  do  not 
provide  possibility  to  decrease  needed  calculations  when 
exists  some  linear  or  non-linear  dependence  between 
given  vectors  X  and  Y.  Proposed  modification  of  NRMG 
algorithm  provides  decreasing  of  needed  calculations 
using  that  NRMG  algorithm  consists  of  sequence  of 
identical  operations  -  multiplications  of  vector  by  Givens 
matrices.  If  some  of  corresponding  components  of  given 
vectors  are  equal,  they  are  not  changed,  which  decreases 
the  number  of  needed  base  operations. 
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